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ABSTRACT 

Our aim in this article is to present some results concerning the inter- 
action of small and large eddies in two dimensional turbulent flows. We show 
that the amplitude of small structures decays exponentially to a small value 
and we infer from this a simplified interaction law of small and large 
eddies. Beside their intrinsic interest for the understanding of the physics 
of turbulence, these results lead to new numerical schemes which will be 
studied in a separate work. 
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INTRODUCTION 


The conventional theory of turbulence in space dimension three implies 
the existence of a length l d which is small in comparison with the macro- 
scopical length lg connected to the geometry, and which is such that the 
eddies of size less than l d are damped by the effect of viscosity and become 
rapidly small in amplitude; the length l d is called the Kolmogorov dissipa- 
tion length [8]. In space dimension two the situation is similar, but l d is 
replaced by the larger length 1 introduced by Kraichnan [9]. It is one 
of our aims in this article to derive directly from the Navier-Stokes equa- 
tions and without any phenomenological consideration a mathematically rigorous 
proof of this property: the exponential decay of the small eddies toward a 

small value. Note however that the cut-off size between small and large 
eddies is much smaller than 1 or even l d , and this is due in part to the 

high level of generality allowed here where singular flows can be considered 
such as those generated by flows in nonsmooth cavities, like the flow in a 
rectangular cavity. A physical discussion on the necessary cut-off length is 
presented hereafter. 

Our approach is the following one: the Navier-Stokes equations of two 

dimensional viscous uncompressible flows are written as 

(0.1) - vAu + (u • V)u — Vw = f in Q x R 

ot + 

(0.2) V • u - 0 in fi x R + 

where u = u(x,t) - {u^ u 2 > is the velocity vector, w = w(x,t) is the 
pressure, f represents volume forces, v > 0 is the kinematic viscosity. As 
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usual (0.1)(0.2) are supplemented by boundary conditions which could be for 
instance 

(0.3a) u = 0 on 9ft 
or 


(0.3b) u • v = 0, vxcurlu = 0 on 9ft, v the unit outward normal on 9ft, 
or 


(0.3c) ft = (0,L 1 )x( 0,L 2 ) and u, w are periodic of period in the 

direction x^, i = 1,2. 

Our emphasis here will be on the space periodic case (0.3c), but the 
other boundary conditions will be considered as well. In all cases (0.1) - 
(0.3) reduces to an abstract evolution equation for u in an appropriate 
Hilbert space H: 

(0.4) 4t" + v Au + B(u) = f. 


The operator A linear, self-adjoint unbounded positive in H with domain 
D(A) C H, is the Stokes operator. Since A - ^ is compact self adjoint, A 
possesses a complete family of eigenvectors w^ which is orthonormal in H 


(0.5) 


Aw j ' x jV 


3 - 1 , 2 , 
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0 < X l £^2»***’ Aj 00 as j -*• ~ 

Of course in the space periodic case (0.3c) the wj are directly related to 
the appropriate sin and cos functions of the Fourier series expansion (see 
[12]). The operator B is a quadratic operator; B(u) = B(u,u), where 
B(» ,• ) is a bilinear compact operator from D(A) into H. 

For fixed m we denote by P ■ P m the projector in H onto the space 
spanned by an< * we wr i te Q = Q m = I - P m » We set 

u=P+q,P*Pu,q = Qu, 

and we show that, after a transient period, and for various norms, p is 
comparable to u, and q is small in comparison with p and u (see Sec. 1). 

We then project equation (0.4) on PH and QH; this yields a coupled 
system of equations for p and qs 

(0.6) || + vAp + PB(p + q) - Pf 

(0.7) iS. + vAq + QB(p + q) = Qf . 

Since q is small in comparison with p one can speculate that B(q,q) = 

B(q) is small in comparison with B(p,q) and B(q,p) and that these quanti- 

« 

ties are small in comparison with B(p,p) = B(p). Also the relaxation time 

for the linear part of (0.7) of the order of (vA ..)”* is much smaller 

nrrl 

than that of (0.6) which is of order (vA^) *. This suggests that an 
acceptable approximation to (0.7) is given by 
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( 0 . 8 ) 


vAq + QB(p) = Qf. 


This leads us to introduce in H the finite dimensional manifold Mg with 
equation 


(0.9) 


q - * n (p) « (vA)" 1 (Qf - QB(p) ) 


p * Pu , q - Qu. 


It is one of our aims to justify this approximation: for large times, i.e., 

after a sufficiently long transient period, the ratio of q to u is of the 
X 1 

order of -r for large m, whereas the distance of q to Mg, (compared 

nH-1 3/2 

to a quantity of the order of u) , is of the order of (■? ) for. 

nH-1 

large m. The proof of this result appears in Sec. 2. Hence, for large time, 
an orbit u(t) ■ p(t) + q(t) corresponding to any solution of (0.4) becomes 
closer to Mg than to the linear space q * 0. We intend in a subsequent 
work to construct a whole family of explicitly defined manifolds provid- 

ing better and better approximations to the orbits as j increases^ ). The 
manifold Mg (as well as the future manifolds Mj ) plays the role of 
approximate inertial manifolds to the two dimensional Navier— Stokes equations 
and constitutes a substitute to them in situations where we cannot prove the 
existence of such manifolds. 

In Sec. 3 we recall and improve significantly a result in [7]: this 

leads us to introduce a Lipschitz manifold £ of finite dimension like Mg» 


< 1 >C. Foias, 0. Manley, and R. Temam, Article in preparation. 


to which all the orbits of (0.4) remain eventually at a distance less than 
exp(-c Hence £ provides a much better approximation than M -. . 

but, on the contrary, the proof of existence is nonconstructive and does not 
provide an explicit expression like (0.9). It provides nevertheless an in- 
teresting complementary aspect. Let us mention also that another type of 
approximate manifold containing all the stationary solutions has been exhibit- 
ed by E. Titi^ 1 ^ . 

This article ends with an Appendix providing a technical but totally new 
method of estimating certain norms of the solutions of an evolution equation 
like (0.4): taking advantage of the analyticity in time of the solutions, we 

estimate the domain of analyticity in the complex time plan and using Cauchy's 
formula, we readily deduce estimates on the derivatives d^u/dt^ from the 
estimates on u in the domain of analyticity; these estimates on the time 
derivatives of u are much sharper than those obtained by real variable 
methods. 

The results presented here were announced in [0]. We intend in a subse- 
quent work to derive approximate manifolds of higher order than Mg and to 
study the three dimensional case. 


1. FAST DECAT OF SMALL EDDIES 

In Secs. 1.1 and 1.2 we briefly recall the functional setting of the 
Navier— Stokes equations and some useful estimates. Then in Sec. 1.3 we derive 
the estimates on the magnitude of the small eddies. 


( 1 ) 


E. Titi, Article in preparation 


1.1 Preliminaries 


As we recalled in the Introduction, the Navier-Stokes equations 
(0.1) (0.2) associated to one of the boundary conditions (0.3) is equivalent to 
an evolution equation 

(i*i) 4^ + vau + * f 

dt 

in an appropriate Hilbert space H. Here f t H, v > 0, A is a linear 

self-adjoint positive operator with domain D(A) C H, and whose inverse A”* 

is compact; we have B(u) = B(u,u) where B(*,«) is a bilinear compact 

operator from D(A) (endowed with the norm |A«|) into H; H is a Hilbert 

2 2 

subspace of L^(ft) . Its norm and scalar product are denoted |»|, (•,•) 

as those of L^(£l)^ or L^(J1); for the details see [11] [12]. 

We recall that for Uq given in H the initial value problem 

(1. 1X1.2): 

(1-2) u(0) = u Q , 

possesses a unique solution u defined for all t > 0 and such that 
(1.3) u e COR^;H)D L 2 (0,T;V), VT > 0; 

here V * D(A^/^) and the norm |A^^» | * on V is equivalent to the 

O 

I- norm of grad u. If Uq e V then 


(1.4) 


u e C( TRj.; V) O L 2 (0,T;D(A) ) , VT > 0 
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In both cases (Uq e H or V), u(.) is analytic in t with values in 
D( A) ; the domain of analyticity of u in the complex plane C t comprises a 
band around 14 . an< ^ is described in more details in the Appendix. 

It is useful here to reproduce some a priori estimates verified by the 
solutions u of (1.1)(1.2). Before that we recall some inequalities (contin- 
uity properties) concerning B (see [7]): for every u, v, w e D(A): 


(1.5) 


( 1 . 6 ) 


I u [ 1/2 lul 1/2 !vl 1/2 lAvl 


1/2 , 1/2 , „ 

u Au II vH 


| B(u ,v) | _< Cj 
| (B(u,v) ,w) | _< c 2 |u| 1/2 luB 1 ^ 2 Ivl [wj 1 ^ 2 . 


where Cj , C 2 like the quantities c^, c^, which will appear subsequently, 
are dimensionless constraints^. Also we recall from [1][3] the inequality 


(1.7) 


L (fl) 


2 - c 3 



*4> e D(A) , 


from which we deduce that 

| B(u ,v) | | (u»V)v| _< 


and using (1.7) 


! |u| . |Vv| 

L (0) 

M M „ 

l (n) 


^ ^ These constants can be absolute constants or they may depend on the shape 
of Jij by this we mean that they are invariant by translation or homothety 
of i]. 
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( 1 . 8 ) 


I ull II vll ( 1 + log — 1- Au -1. -) 

X ^ II u II 1 

I *3/2 ,2 1/2 

|u| |Av| (1 + log -L- — yi-) 


2 1/2 


B(u,v) | £ c 4 


1.2 Behavior of Small Eddies 

As mentioned in the Introduction we fix an integer m e N and denote 
by P = P m the projector in H onto the space spanned by the first m 
eigenvectors of A, Wj,»»»,w ; we set also Q = Q m = I - P m , and for the 
sake of simplicity 


(1.9) 



A » X 


m+1 * 


We write p - Pu , q * Qu; p represents a superposition of "large eddies" of 

size larger than x”**^ 2 , and q represents "small eddies" of size smaller 

m 

than X~i{ 2 . By projecting (1.1) on PH and QH we find since PA * AP 
m+1 

and QA = AQ: 


(1.10) + vAp + PB(p + q) = Pf 

(1.11) + vAq + QB(p + q) “ Qf . 

We take the scalar product of (1.10) with q in H: 

(1.12) |q | 2 + vU qll 2 = (Qf,q) - (B(p + q) , q) . 


Thanks to the orthogonality property 



(1.13) 


(B(<|> = 0, 


*♦, + e V, 


the right hand side of (1.11) reduces to 


(Qf,q) - (B(p,p),q) - (B(q,p),q). 


Using (1.6) and Schwarz inequality we majorize it by 

■ .2 1/2 

I Qf | |q | + c.apir |q| (1 + log - L - P -L ) + c 9 | q | iqi i P n 

x 1 ipii 

(since II pll II ull ) 

2 1/2 

_< | Qf | | q | + c. » pP 2 |q| (1 + log -L A - P -1__ ) + c„A" 1/2 flql| 2 Hull . 

X 1 lp» 

We denote now by M 0 (resp. Mj, H 2 ) a bound of |u| (resp. | U |, |Au|), on 
the interval of time I = (t^ ~) under consideration 


(1.14) M- - Sup|u(s) | , M » Suplu(s) I , M = Sup|Au(s)|; 

sel sel sel 


we observe that 

I Ap | 2 _< X m » pi 2 = XI p* 2 


and set 


(1.15) 


L = (l + log 


m+1 > 


We obtain 
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(1.16) ^ | q | 2 + (2v - c 2 A 1 / 2 Mj) «q « 2 < |Qf| |q| + c 4 M 2 L 1 / 2 |q| . 

-1/2 

Hence, assuming that c 2 A <, v, i.e., 

2 c 0 M. 2 

(1.16) yields 

| q | 2 + I' q" 2 < A _1/2 ( | Qf | + (AM 2 L 1/2 ) «qll 
<7 »q « 2 C iQf | 2 + (c 2 mJl) 

(1.18) I q 1 2 + v n qn 2 < ^ ( |Q f 1 2 + c 2 m^l) 

(1.19) lq|2 + vA ' q i 2 i^A ( l Qf|2 + C 4 M 1 L) - 


We infer easily from (1.19) that for t t ^ , t e I: 

(1.20) ( q( t) f 2 < f q( t j ) 1 2 esp(-vA(t - tj)) + C f Qf | 2 + c 2 MjL). 

v A 

Before interpreting this inequality, we derive a similar inequality for the 
(Hi) V norm. Taking the scalar product of (1.11) with Aq in H we find 


7 11 ql1 2 + v | A-q | 2 = (Qf ,Aq) - (B(p + q) , Aq) . 


We expand and use Schwarz inequality and (1.6) - (1.8) to majorize the right 
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hand side of this equation by 


| Qf | | Aq | + c-ll pllL 1/2 | Aq | (II pll + llqll) 


+ c, |q | 1/2 | Aq | 3/2 ( II pll + II qn) 


< (with Young's inequality) 


<? |Aq| 2 |Qf| 2 

V 


Thus , 


( 1 . 21 ) 


( 1 . 22 ) 


4 2 4 

Ij- iql 2 + V |Aq| 2 < c' 3 (i |Qf| 2 + -1- + 


4 2 4 

. , 9 M^L 

^ llqll 2 + vMlqll 2 < c'(i | Qf | 2 + — + -3*} 


and we conclude that 


It q( t )ll 2 < II q(t 1 )ll 2 exp(-vA(t - t^)) 


(1.23) 


m?h! 


+ ^ ( i | Qf | 2 + 3. L + ^l) 

vA v v 1 1 v 3 J 


. 9 0 9 

In (1.20) and (1.23) we can bound IqU^ 2 and B q(t x )I1 by M Q and Mj^ 
respectively. Then after a time depending only on Mg (or M^), v and 
A = A j, the term involving t becomes negligible and we obtain 
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(1.24) 


|q(t) | 2 <-%-? (|Qf | 2 + c 2 m*l), 
v A 

# q(t)ll 2 < -^ ( I Qf | 2 + M^L + -^i-) 

v A v 


for t large. Alternatively, denoting by k , k', some quantities 

which depend only on the data v, f, ft, and Mq, Mj, M 2 » we rewrite (1.24) 
as 


(1.25) 


|q( t) | 2 < k L6 2 , Hq(t)l 2 _< < L6 for t large, 


. X i X 1 T . . i X m+1 

5*7— = 7 , L - 1 + log 


A X 


nH-l 


X 1 * 


Using also the results in the Appendix we conclude the following 


Theorem 1.1: We assume that m is sufficiently large so that (1.17) 

holds. Then for any orbit of (1.1), after a time t* which depends only on 
the initial value u(0) = uq, the small eddies component of u,q = Q m u > is 
small in the following sense 

| q( t) | <k q L 1/2 S, »q(t)B < Kl L 1/2 6 1/2 

(1.26) |q'(t) < k'L 1/2 6 

| Aq( t) | £< 2 l1 ^ 2 ’ t 2 **• 

The first two inequalities in (1.26) follow from (1.25); the third one follows 
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from (1.25) and the analog of (A. 15) for q ^D. The fourth inequality is 


obtained by writing 


vAq = Qf - q' - QB(p + q) 


|Aq| <i | Qf | +1 |q'| +1 |QB(p + q) 


and utilizing (1.5), (1.6), (1.8). 

In Sec. 1.3 hereafter we intend to provide a more explicit form of the 
constants k in the case of space periodic flows. 


1.3 The Space Periodic Case 

We first review the well-known a priori estimates on the solutions of 
(1.1). This will yield more explicit expressions for Mq, Mj, M 2 » 

We take the scalar product of (1.1) with u in H; using the ortho- 
gonality property (1.13) we obtain 


1 d . ,2 


u| + vdull = (f ,u) _< | f | |u 


< x x 1/2 | f | null 


(1.26) 


<T» u » 2+ 2^1 f 1 


d | u |2 +V | un 2 < 1 | f 1 2 


Note that q is analytic in the same region of the complex plan as u. 
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Theorem 2.1: For t sufficiently large, t > t^, any orbit of (1.1) 

1/2 

remains at a distance in H of PH of the order of kL 6 and at a 

— m ■■■■■■— - 

distance in H of Mq of the order of K In the norm of V, the 

1/2 1/2 

corresponding distances are of order <6 L and «L6 ; the constants 

k depend on the data v, |f|, and t* depends on these quantities 

and on Rq, when |u(0)| £ Rq. 

3. A NONCONSTRUCTIVE RESULT 

Our aim in this last section is to exhibit a manifold £ which is 

Lipschitz, has finite dimension, and captures the solutions of (1.1) in a much 
narrower neighborhood than Mq does. However, the existence of £ is. 

proved in a nonconstructive way, by opposition with the very simple and ex- 
plicit equation (2.2) available for Mq. Secs. 3.1 and 3.2 provide prelimi- 
nary results and Sec. 3.3 contains the main one. 

3.1 Quotient of Norms 

We consider two solutions u, v of (1.1) and set w » u - v: 

(3.1) ■JF + v Au + ®( u ) “ u(0) * u 0 > 

(3.2) + v Av + B ( v ) * v(0) * v Q , 

(3.3) + v Aw + B(u,w) + B(w,v) ■ 0. 

at 

2 2 

Let a denote the quotient of norms Iwl /|w| ; then 
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da 2((w',w)) 2llwll , ~ \ 2 . , s 

•jr - 5 7— (w ,w) = T (w , Aw - aw) 

dt |»| 2 Nr Iwl 2 


j (\)Au + B(u,w) + 3 (w,v), Aw - aw). 

Iwl 


2 

Since (Aw, Aw - aw) = Aw - aw , we conclude, using ( 1 . 5 ), that 


da , 2v 


+ — | Aw - aw 


dt 1 1 2 

w 


it (B(u,w) + B(w,v) , Aw - aw) 

Iwl 2 


£ 2 2 |Aw - aw| ( |B(u,w) I + |B(w,v) | ) 


w 


2 c , 


< ^ I Aw ” &w| 2 ( |u| 1/2 I Au 1 1 ^ 2 H wll + I w| 1/2 |Aw|^ 2 IIvB) 

|w| 

2 c 2 

V 2 I Aw - aw | 2 + ( | u | | Au | + II vll | Av | X ^ ^ 2 )a . 

|w| V 


Hence 


( 3 . 4 ) 


da 

dt 


Aw - aw | < pa 


w 


where 

P - — j-tj lul I Au | . 
vX l 

By integration of the differential inequality a" < p a, we find that for 
tj<t<T<tj + T 


( 3 . 5 ) 


P + P • 
U V 
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/o II w(t )ll 2 , Hw(t)U 2 /f T f \ j \ 

(3.6) — - ■ - • o ' < — _ exp(J p(s)ds). 

|w(x)| 2 " |w(t)| 2 t 

Now we estimate the integral of p in terms of the data; as in (1.16) 
we assume that on the interval of time under consideration 

(3.7) llu(t)ll<M 1 , lv(t)l£M 1 . 

With an appropriate value of Mj (3.7) will be valid on some finite interval 
of time [0,T], or on some interval of time (t^, 00 ), once the orbits have 

entered the absorbing set. 

We have 


t 2c. T 

/ P u ds £ — i7J I llu!l l Au |ds 
t vX j t 

<4„(,. t) W (/ T |Au| 2 d8, l/2 . 
vxj' t 

An estimate on Au is obtained by taking the scalar product of (3.1) with 
Au in H: 

Hull 2 + 2v | Au | 2 =* -2(B(u),Au) - 2(f,Au) 

<. 2 1 B(u) | | Au | + 2 | f | | Au | 

_< (with (1.5)) 

<_ 2c ^ |u| 1/2 Hull | Au | 1 ^ 2 + 2 | f | | Au | 
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_<v|au | 2 + -| |u | 2 a ufl 4 + ^ I f | 2 
v 


(3.8) 


^lui 2 +v|Au| 2 <i |f | 2 +-^- M r 

v A, 


Thus , 


(3.9) 


/ |Au| 2 ds < ^ (|f| 2 + -J- Mj) 


V X . 


t ! +T m 2 0 

1 o M. « r |2 £ 

/ |Au| 2 ds < -i + lLL + -J-_ M* 

t. v v X, 


and 


(3.10) 


t 
c 
vX 


T , 

/ p u ds < Tj- (-r - t) 


1 / 2 .. 


M^" 9 1/2 

' -1 v v *1 


Since the estimates on v and p are the same, we have 

v 


(3.11) 


/ p ds < (t - t) 1 ^ 2 K_. 
t J 


3.2 The squeezing property 

The squeezing property is an important property of the solutions of the 
Navier-Stokes equations which has been introduced in [6]. A stronger form of 
it, called the strong squeezing property or the cone property was proven in 
[4] for some other, more strongly dissipative equations. For the two dimen- 
sional Navier-Stokes equations, we derive here a form of the squeezing proper- 
ty sharper than in [6]. 
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We take the scalar product of (3.3) with w in H and thanks to (1.13), 


(1.16) we find 


|w| 2 + 2v II wll ^ = -2b(w, v, w) 


< 2c 2 |w| II wll llvll 


< V II wll 2 + — I w I 2 II vll 2 
— V 1 1 


O c o O O 

< vll wll + — M. I w I 

— V 1 1 1 


(3.12) 


2 c 

d | 1 2 , II wll l 2 Jn, .2 „ n 

_ | w | + (v — 7 - — Mj)|w| < 0. 


consider tp, t, 0 < t < tp < T and write, using (3.6)(3.11) 


(3.13) 


Y 0 = 


— )) " 2 < exp(K (t - t) 1/2 ) . 

|w(t 0 )| 2 " 3 ° |w(t) 1 2 


Thus, 


( 3 .!4) 1- |w| 2 + (vy Q exp(-< 3 t* /2 ) - ^ M^) |w| 2 _< 0 


1 / 2 , l 2 


and by integration 


(3.15) |w(t Q )| Z < |w(0) exp(-VY 0 t Q exp(-K 3 tJ /2 ) + ^- M 2 t Q ) 


Now if |Q m w(t 0 ) | > |P w(t 0 )| , we write 
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and 


(3.16) 


|P m w(t 0 )|2 + " Q m w(t 0 )l12 

° l p m w(t o ) l 2+ IV^o^ 2 


2 l Q m w(t O ) 


Vi 


! w(t O ) I 2 — l w(0) | 2ex P (_vX nr fl < 5 t O + lc 4 t O ) 


K , = — 

4 v 


2 2 1 , 1 / 2 , 

~ M l’ K 5 = I exp ^ -K 3 t 0 


Of course the interval (0,t 0 ) can be replaced by any interval (tj , tj + t Q ) 
on which the bound (3.7) is valid. 

In conclusion (this is the squeezing property ), whenever (3.7) is valid 
on some interval (tj, tj + tQ)» then w a u — v satisfies one of the 
following conditions: 


(3.17a) 



or 

(3.17b) |w(t Q + t x )| 2 < | w(tj) | 2 exp(-vA mfl < 5 t () + 

Since < 5 are independent of m, the exponential term in (3.17b) can be 

made arbitrarily small by choosing m sufficiently large; we will take advan- 
tage of this remark in Sec. 3.3. 

Of slightly more explicit form of can be derived by using the 

Grashof number G - |f|/v 2 X 1 and the Reynolds type number = M^vA^ 2 . 


We find (t ■ t^): 
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K 3 - c^Cvx/'X ♦ VX,G 2 + t 0 vX lR ‘) 1/2 

(3.18) = c 2 R n^ vX l^ 

k 5 “ 1 e *P(<- c ?„<' , Vo ,I/2 ° t n + V X l c2 + V X l R n )1/2 )- 

In the space periodic case we have seen that, for large times, we can take Mj 
= (2|f|G)l/2. Then = FL G and the above quantities become 

< 3 = c'(vX 1 ) 1 / 2 (G 4 + t^G 4 + t 0 vX lG 8 ) 1/2 

(3.19) < 4 = 2c 2 (vX 1 )G 2 

k 5 = T + ^^l 0 * + tQvXjG 8 ) 1 ^ 2 ). 

3«3 The Approximate Manifold 

We denote by S(t), t > 0 the operator in H: + u(t), where u(*) 

is the unique solution of (1.1) satisfying u(0) = ug* The operators S(t), 
t 0, form a semigroup in H. 

The squeezing property tells us that if u(*)» v(0) are two solutions 
of (1*1)) lying in the ball (<J> e v, li< M < M^}, for 0 t _< T, then at 

each time t e [0,T] and for every m e we have either 

|Q n (S(t)u 0 - S(t)v 0 )| < |P m (S(t)u 0 - S(t)v 0 )| 
or 
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|S(t)u 0 - S(t)v Q | < |u 0 - v Q | exp \ (-v^lVo + < 4 t 0 ) 
as above. 

Now we choose tQ e [0,T], m e H and consider a subset I " I (m) of 

S(t Q ) { u Q e V, llu 0 ll < Mj} 

which is maximal under the property 

(3.20) |Q (u - v)| < |P (u - v) | . 

m — m 

By this we mean that if u e ^ (m) then 

{v e V, v satisfies (3.20)} £ (m). 

The existence of such a maximal set Is easy. 

We then apply the squeezing property: whenever II u ( s ) II , we see 

that S(to)u(s) = u(t 0 + s) either belongs to £ (m) , i.e., 

|Q m (S(t 0 ) u(s) - S(t 0 H)| < |P m (S(t 0 )u(s) - S ( 1 0 )4> ) | , 

for some 4> e V such that H<t>U < and S ( tQ )4> e £ (m) or, if not, 

then for every such <f> 


|S(t Q )u(s) - S(t Q )(J> | 2 < | u(s) - <))| 2 ex P( -vX nH . 1 lc 5 t o + K 4 t 0^ 
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< _ expC-vX^jKgtj, + < 4 t 0 ). 


In all cases the distance of S(tg)u(s) to ][ (m) is bounded by 


■ M i fc 0 

172 exp (— (ic 4 ~ vX m +l K 5 ) ) 


We can choose t^ = (vX^) and the bound becomes 


provided that 


2M i K c X 

1 , 5 m+l N 

7172 ex P ( - T T-) 

" 1 1 


(3.21) 


X ro+1 > 2k 4 

Xj — k 5 v ^i 


By translation in time (t ■+■ t — t^), we conclude that once the orbit u has 
entered the absorbing set { n 4> II M } , which happens for t > t* = t+(R n ) 

(for | u( 0 ) | < Rq), the distance of S(t)ug to £ (m) is bounded by a 
given quantity E, 


(3.22) 


dist H (S(t)u 0 , l (m)) < E 


provided t>t. +(vX.) , and 


i • e • y 


, *5 Vl. . _ 

exp ( " T —■ 1 < E ' 


4 - _ 

-T > - — log E. 

A 1 " k 5 


(3.23) 
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By definition the set £ (m) enjoys the property that 

|Q (u - v) | < |P (u - v) | , Vu, v e M(m). 
m — m 

Hence, £ (m) is the graph of a Lipschitz function 

<p: P l (m) ♦ OH 

m m 

|f(P u) - f(Pv)| < |Pu - P v|, V P u, P v P m l (m). 
m m — ra m mm m 

By the Kirszbaum extension Theorem [16] ¥ can be extended as a Lipschitz 

function (with the same constant) form P m H into Q m H, that we still denote, 
by Now ¥ is defined from P m H into Q m H> and its graph is a 

Lipschitz manifold above all of P m H. 

In conclusion we have proved the following theorem 
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APPENDIX 


Estimates in the Complex Time Plane 

It was proved in [6] (see also [12]) that the solutions to the Navier- 
Stokes equations are analytic in time; we want to show how one can then use 
Cauchy^s formula to get a priori estimates on the time derivatives of the so- 
lutions# The main point in the proof is to determine the width of the band of 
analyticity of the solution around the real axis Il|_; this will follow as 

in [6,12] from a priori estimates on the solution in the complex plan. 

Id 

The complex time is denoted £ = se ; m F, 1D(A) are the 

complexified spaces of H, V, D(A); A, B are extended as linear and bilinear 
operators respectively from 1D(A) into TH : 

(A.l) A(u^ + iu2) = Auj + iAu2» 

(A. 2) B(uj + iU2» Vj + iv2) = B(uj,Vj) - B(u2»V2) 

+ i[B(u2»vj) + B(uj,V2)1 

Vu = + iu^, v = + iv£ e H(A). The Navier-Stokes equation (1.1) 

becomes (u = u(c)): 

(A. 3) + vAu + B(u) = f 


(A. 4) 


u(0) = uq 
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i oo A 

Assuming that uq V(or ^E/), u|j^eL (I^; and as in (1.14), we denote 

by Mq, M lf the supremum of |u(t) | and II u( t ) II » t e 3R,|_. We take the 

scalar product in H of (A. 3) with Au; we multiply the resulting equation 
i0 

by e and take its real part. This yields 

II u(se^ )ll ^ + vcos0 |Au(se^)|2 = 

(A. 5) = -Re e* 9 (B(u),Au) - Re e* 9 (f,Au) 

_< | (B(u) ,Au) | + |f | |Au| . 


We expand by bilinearity (using (A. 2)) and bound the resulting expressions 
with the help of (1.8): 


9 I a I 2 1/2 

|(B(u),Au)| < cllull (l + log I U l a) | Au | . 


A j Hull 


Also 


I Au | < 


VCOS0 


- 2 


|Au | 2 + 


2vcos0 


Hence (with u 


u(se i9 ) ) : 


(A. 6) 



+ vcos0 Au 


< 


< ■ + e 

— vcos0 c 5 


2 

Hull | Au | (l + log 


| Au 
X j B till 


2 1/2 

• 


We write 



and consider the function 
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XjVCOS0 2 \/2 21/2 

z + <|> ( z) z + c,. llullX z( 1 + logz ) , 


By elementary computations 


( 1 ) 


(A. 7) 


2 2 2 2 
Oull 4c^«uir 

♦ (z) i 2uco ' s 0 log (: 2 for 2 l 1 * 


2 2 
XjV cos 0 


and (A. 6) yields 


d . 2 , vcos6 | A 1 2 , 
__ „ uU + __ | Au | < 


(A. 8) 


2 2 
4c_HuU 


I £ 1 4 c. . 4c,Huy 

+ ( lp g ~ g Z~ ) • 


Setting 


X^v cos 0 

(1 + 4c 2 ) . ^ 

y(s) * ^ (|f| + #u(se °)( ) we infer from (A. 8) that 

XjV cos 0 


lii c i ^c 083 ® y 2 lo S y , 


where c' is an appropriate nondimensional constant. As long as 
y(s) 2y^ * 2y(0), we have 


( i ) 2 2 2 2 

Looking for the maximum of -az +6(1+ log z ), we find 


0 2 ( 1 + log z 2 ) < a 2 z 2 + 0 2 log ~ 


z0(l + log z 2 ) 1/2 < az 2 + Bz(log 2j) 1/2 

a 

.2 ,2 


i 2az 2 + \ (log ^). 


We then choose a * 


X i vcose 1/2 

-^“4 , 8 - c 5 luixj /2 . 
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y' < c'XjVcosQy log(2y Q ) 


— 1 - Cj X 1 vcos61og(2y 0 )s ’ 


and this is indeed 2y^ as long as s _< T^: 


T = X - 

* 2cjX 1 vcos9y 0 log(2y 0 ) 


For #UqII _< Mp we replace T* by 


(A. 9) T * (1 V " 


ZcJXjVcosef-Ej- + _ T ‘)l„g 2 (-V + 2 * » ) 


COS 0 XjV COS 0 


COS 0 XjV COS 0 


(A. 10) 


u(se 10 )n 2 < 2 ( | f | + nu 0 n^) < 2(|f| + m*) 


0 < s < 


3cos0 


M. r n 

2c'X 1 v(G + j) + log 2( j- + 2 jp) 


COS 0 XjV cos 0 


and in particular for 
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CA.ll) 


0 < s < 


3cos0 

m" 2 " 


M 1 M T 

2cJX 1 v(G + — y + log 4(G + ? ) 


X 1 V 


X 1 V 


when cos^Q * 

Following the method developed in [6] we conclude that the solution u 
of (A. 3) (or (l.a)) is analytic in the region 


(A. 12) 


i0 

A(Uq) = {se 1 , s _< a cos 0, cos 0 } 


a = 


M, 


Mt 


2c'X lV (G + — L-) + log 4(G + — ^ •) 
which comprises the regions 


X ! v 


X 1 V 


| Im C | < Re ( , 0 < Re C _< j 


and 


(A. 13) 


1 1m £ £ » Re ? • 


At any point t e 1^, t ^a, we can apply Cauchy's formula to the circle 


T centered at 

t of radius < 

x/4: 

(A. 14) 

d k u( t) 
, k 

3 15 

It 


dt 

r 

Thus , 



(A. 15) 

Sup 

d k u(t) 

If 


t>a 

dt k 


u(c) 


r (t-c) 


k+1 


dC 


i T k! M 0 

a 
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(A. 16) 


1c 1c 

Sup ll^-£ (t)ll < k! Mj, 
t>a dt a 


Explicit values of Mq and were derived in (1.36) for the two 

dimensional space periodic case: = (2|f|G)^^(t This yields 

(assuming G 2 1)* 


a * 


2c'X 1 v(G + 2G 2 )log 4(G + 2G 2 ) 


( A. 17 ) 


a > 


XjVG logG 


and we deduce from (A. 15), (A. 16) that for t sufficiently large 


( 1 ) 


(A. 18) 


1 


|AoO| < c If 1/2 

dt 

£ 

dt 


j fyy- (|f|A 1 ) k/2 (G 2 logG) 1 


< c|f| 1/2 (|f|X 1 ) k/2 (G 2 logG) k . 


In particular (k = 1): 


(A. 19 ) 


l^r 1 ! < c|f| tW 


| d “^. t> l < c|f|x[ /2 G 2 logG, t > T,. 


This produces an interesting bound on |A u (t)| for t large: 


^^This means as in Theorem 1.1 and elsewhere t > T*(Rq,v , Xj , | f | ) » for 

I u 0 1 - R 0 * 
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(A.2 0) 


vAu = f - B(u) - u ' 

l Au l < 7 1^1 + ^ M 1 / 2 .u. |Au | 1/2 +1 ju'j 

2 

|Au| <1 f f | + ~L f u f null 2 +1 Ju'| 

v v 

< c( I f |x 2 ) ly ^ 2 (G 1>/2 + G + G 5 // 2 logG) 

| Au(t) | < c(|f Ixp^V^logG, for t > T*. 


-40- 


REFERENCES 

[1] H. BrSzis and T. Gallouet, "Nonlinear Schroedinger evolution equa- 
tion," Nonlinear Analysis Theory Methods and Applications , Vol. 4, 
1980, p. 677. 

[2] C. Foias, 0. Manley, and R. Temam, "Sur l'interaction des petits et 
grands tourbillons dans des Scoulements turbulents," C. R. Ac. Sc. 
Paris , 1987. 

[3] C. Foias, 0. Manley, R. Temam, and Y. Treve, "Asymptotic analysis of 
the Navier-Stokes equations," Physica 6D , 1983, pp. 157-188. 

[4] C. Foias, B. Nicolaenko, G. Sell, and R. Temam, "VarietSs inertielles 
pour V equation de Kuramoto-Sivashinsky ," C. R. Ac. Sc. Paris , 301, 
S§rie I, 1985, pp. 285-288 and "Inertial Manifolds for the Kuramoto- 
Sivashinsky equations and an estimate of their lowest dimension," J. 
Math. Pure Appl. , 1988. 

[5] C. Foias and G. Prodi, "Sur le compartement global des solutions non 
stationnaires des equations de Navier*-Stokes en dimension 2," Rend . 
Sem. Mat. Padova , Vol. 39, 1967, pp. 1-34. 

[6] C. Foias and R. Temam, "Some analytic and geometric properties of the 
solutions of the Navier-Stokes equations," J. Math. Pure Appl* , Vol# 
58, 1979, pp. 339-368. 



-41- 


[7] C. Foias and R. Temam, "Finite parameter approximative structures of 

actual flows," in Nonlinear Problems: Present and Future , A. R. 

Bishop, D. K. Campbell, B. Nicolaeuko (eds.), North Holland, 
Amsterdam, 1982, 

[8] A. N. Kolmogorov, C, R. Ac. Sc. URSS , Vol. 30, 1941, p. 301; Vol. 31, 
1941, p. 538; Vol. 32, 1941, p. 16. 

[9] R. H. Kraichnan, "Inertial ranges in two dimensional turbulence," 

Phys. Fluids, Vol. 10, 1967, pp. 1417-1423. 

[10] G. MStivier, "Valeurs propres d'opirateurs dSfinis sur la restriction 
de systemes variationnels a des sous-espaces ," J. Math. Pure Appl. , 
Vol. 57, 1978, pp. 133-156. 

[11] R. Temam, Navier-Stokes Equations , 3rd Revised Ed., North Holland, 

Amsterdam, 1984. 

[12] R. Temam, Navier-Stokes Equations and Nonlinear Functional Analysis , 

NSF/CBMS Regional Conferences Series in Appl. Math., SIAM, 

Philadelphia, 1983. 

[13] R. Temam, Infinite Dimensional Dynamical Systems in Mechanics and 

Physics, Springer-Verlag, 1988. 


-42- 


[14] J. H. Wells and L. R. Williams, Imbeddings and Extensions in Analysis , 
Springer— Verlag, Heidelberg, New York, 1975. 


(VASA 

National Aeronautics and 
Spade Administration 


Report Documentation Page 


2. Government Accession No. 


1. Report No. 

NASA CR-178370 
ICASE Report No. 87-60 


4. Title and Subtitle 

ON THE INTERACTION OF SMALL AND LARGE EDDIES 
IN TWO DIMENSIONAL TURBULENT FLOWS 


3. Recipient's Catalog No. 


5. Report Date 

September 1987 


6. Performing Organization Code 


7. Author(s) 


C. Foias, 0. Manley, R. Temam 


8. Performing Organization Report No. 

87-60 


10. Work Unit No. 

505-90-21-01 


11. Contract or Grant No. 


13. Type of Report and Period Covered 


14. Sponsoring Agency Code 


9. Performing Organization Name and Address 

Institute for Computer Applications in Science 
and Engineering 

Mail Stop 132C, NASA Langley Research Center 
Hampton, VA 23665-5225 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, VA 23665-5225 


15. Supplementary Notes 

Langley Technical Monitor: 
Richard W. Barnwell 

Final Report 


16. Abstract 

Our aim in this article is to present some results concerning the inter- 
action of small and large eddies in two dimensional turbulent flows. We show 
that the amplitude of small structures decays exponentially to a small value and 
we infer from this a simplified interaction law of small and large eddies. 
Beside their intrinsic interest for the understanding of the physics of turbu- 
lence, these results lead to new numerical schemes which will be studied in a 
separate work. 


17. Key Words (Suggested by Author(s)) 
fluid mechanics, turbulence, eddies 
simulation 


18. Distribution Statement 

02 - Aerodynamics 
34 - Fluid Mechanics and Heat 
Transfer 


19. Security Classif. (of this report) 

Unclassified 


NASA FORM 1626 OCT 86 




20. Security Classif. (of this page) 
Unclassified 


21. No. of pages 

44 



NASA-Langley, 1987 















